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' Abstract 

o 

' The leading twist light-cone distributions for transversely polarized p-, p - and 61-mesons are 

re-analyzed in the framework of QCD sum rules with nonlocal condensates. Using different kinds 
of sum rules to obtain reliable predictions, we estimate the 2-, 4-, 6- and 8-th moments for 
transversely polarized p- and p'-meson distributions and re-estimate tensor couplings /J , b± . We 
stress that the results of standard sum rules also support our estimation of the second moment 
C*~) . of the transversely polarized p-meson distribution. New models for light-cone distributions of 

these mesons are constructed. Phenomenological consequences from these distributions are briefly 
discussed. Our results are compared with those found by Ball and Braun (1996), and the latter 
\^ ' is shown to be incomplete. 

o ' 
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^ 1 Introduction 



In this paper, we complete our investigation of the leading twist light-cone distribution amplitudes 
(DAs) for lightest transversely polarized mesons with quantum numbers J PC = 1 (p±, p' ± ), l" 1 

in the framework of QCD sum rules (SRs) with nonlocal condensates (NLC). These DAs are 
important ingredients of the "factorization" formalism |Q] for any hard exclusive reactions involving 
p-mesons. For this reason, the DAs have been attractive for theorists for a long time: the main points 
are presented in |@, @] , a detailed revised version of the standard approach is in j|] , and a generalization 
to the next twists is in [||. The leading twist DA ip T p , bi (x, p?) parameterizes the matrix elements of 
the tensor current with transversely polarized p(770)- and /?'(1465)-mesons (J = 1 ) 

(0 I u(*)<V d (°) I P±(P> A )>L 2=0 = */p1 MP, VP* - e v (p, A)p M ) J dx e ix ^ ^J ± (x,p 2 ) + . . . , (1) 
and the 6i(1235)-meson 

(0 I u(z)a^d(0) I b+(p, A)) = ftfy„ aP s a (p, A)/ /' dx e ix ^ (x, p 2 ) + . . . (2) 

z — " J 



(here dots represent higher-twist contributions, explicitly defined in Appendix A, see Eqs.( [A.9| )-( A.10| ) 
and ref.||). In the above definitions, p v and e^ip, A) are the momentum and the polarization vector 
of a meson, respectively, and p? is normalization point. 
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In the framework of the standard approach, one should restrict oneself to an estimate of the second 
moment (£ 2 ) of the DA to restore its shapeQ. In other words, the variety of different DA shapes is 

r o /o 

reduced to the 1-parameter family of "admissible" DAs: ip(x;a2) = 6x(l — x) 1 + a>2C 2 (2x — 1) . 
This family includes both the asymptotic DA (02 = 0) and Chernyak-Zhitnitsky model || for the 

7r| C2 

pion DA (aj 1 = -2/3). For the pion case, one can think it is rather enough: most of debates (see 
[§, ||, [H], 12] and refs. therein) about the shape of this DA are concerned just with the value 

7T I CZ 

of coefficient a<i - is it close to or to a 2 ? In our opinion, advocated since 1986 @, the shape 
of the pion DA is not far from the asymptotic one || ^, [l3|, 14, ||]. Only recently, researchers have 



tried to extract the next Gegenbauer coefficient jl2| and other parameters of the pion DA [15] from 
experimental data. But, in general, there is no principle to exclude a more rich structure for a hadron 
DA. In this case, the standard approach is definitely out of its applicability range, and one should 
use more refined techniques, e.g., the QCD SRs with NLC. 

This work was started in || where the "mixed parity" NLC SR for DAs of p- and 61-mesons, the 
particles possessing different P-parity, was analyzed. We concluded that, to obtain a reliable result, 
one should reduce model uncertainties due to the nonlocal gluon contribution. Separate SRs for each 
P-parity channel should be preferable for this purpose, and here we construct these "pure parity" 
SRs for corresponding DAs. The SR of this type possesses a low sensitivity to the gluon model but 
involves contributions from higher twists^. To construct a refined "pure parity" SR for twist 2 DA, 
one must resolve the corresponding system of equations (see Appendix A). We realize this solution 
using the duality transformation, introduced in our previous work Eg ]. The negative parity NLC 
SR for the transversely polarized p-, //-mesons works rather well and allows us to estimate the 2-nd, 
4-th, 6-th, and 8-th moments of the leading twist DAs. The positive parity SR for the transversely 
polarized 61-meson can provide only the value of the 61-meson tensor coupling, ff. We suggest the 
models for these DAs and check their self-consistency, based upon both "pure" and "mixed" NLC 
SR. The DA shape ipj ± (x) differs noticeably from the known one. Finally, we inspect how these 
models can influence the B — * pev decay form factors. 

The approach has been grounded in [DJ , the calculation technique is the same as in |7], || ; 
therefore, the corresponding details are omitted below. Some important features of the NLC SRs 
approach would be briefly recalled. The original tools of NLC SR are nonlocal objects like Mj(z 2 ) = 
(q(0)E{0, z)q(z)} or M$(z) = {q(0)^E(0, z)q(z)), rather than constant quantities of (g(O)g(O))- 
type. Note that, in deriving sum rules, one can always make a Wick rotation and treat all the 
coordinates as Euclidean, z 2 = —z\^ ucl < 0. NLC Ms(z 2 ) can be expanded in the Taylor series over 
the standard (local) condensates, (q(0)q(0)), (<?(0)V 2 (/(0)), and over "higher dimensions" (see details 
of the expansion of different NLCs in ]p^| ), 

M S (z 2 ) = (g(O)g(O)) - ^(q(0)V 2 q(0)) + . . . . (3) 

So, one can return to the standard SR by truncating this series. But, in virtue of the cut off, one 
loses an important physical property of nonperturbative vacuum - the possibility of vacuum quarks 
(gluons) to flow through vacuum with a nonzero momentum k q ^ 7^ 0. The parameter (k 2 ), fixing 
the average virtuality of vacuum quarks, can be interpreted as a measure of condensate "nonlocality" 

\2 



(g(O)g(O)) " 2(g(0)?(0)) l^hiral hmitj 



1 We should note in this respect that the standard approach could not provide a reliable estimate even for the second 
moment of DA, see [0, Q, §, § 
2 as was noted in H. 

3 Here E(0, z) — Pexp(i J dt^A^^Ta) is the Schwinger phase factor required for gauge invariance. 
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The A 2 was estimated from the mixed condensate of dimension 5, A 2 ~ 0.4 — 0.5 GeV 2 [19, 20|. It 
is important that its value is of an order of the characteristic hadronic scale, A 2 ~ m 2 ~ 0.6 GeV 2 , 
therefore the nonlocality effect can be large, and it should be taken into account in QCD SR. Really, 
the second term in the expansion (^|) of M$(z 2 ) that is inverse of the first one in sign becomes of 
an order of the first term at \z 2 \ ~ 1/m 2 due to the estimate |A 2 z 2 | ~ 1. Moreover, we should take 

into account the whole set of ^A 2 z 2 ^ -type corrections, appearing in the Taylor expansion. These 
corrections just form the decay rate of the NLC (Ms(z 2 )) in the main. The sensitivity to this rate is 
crucial for the DA moment SR: it leads to much softer behavior of DA near the end points x = 0, 1 
and allows one to extend QCD SR to higher moments (£ N ) = Jq (p(x)(2x — l) N dx, as it was shown 
in §. 

Since neither QCD vacuum theory exists yet, nor higher dimension condensates are estimated, it 
is clear that merely the models of NLC can be suggested (Appendix B). Here we apply the simplest 
ansatz to NLC |8| that takes into account only the main effect (k 2 ) = A 2 ^ and fixes a length 
of the quark-gluon correlations in QCD vacuum A = l/\ q ~ 0.8(Fm) ||, M. This suggestion leads 
to the simple Gaussian decay for Ms(z 2 ), while the coordinate behavior of other NLCs looks more 
complicated. Certainly, the single scale of decay for all types of NLC (see Appendix B) looks as 
a crude model. But, the model can be rather crude if one deals with SRs only for the first few 
moments (£ N ), because for these integral characteristics the details of NLC behavior appears to be 



not very important (see discussion in sect. 5). An alternative case is provided by a special SR []13|, |14] 
constructed directly for the shape of DAs. 

In nowadays, the lattice calculations of NLC provide an inspiring knowledge pl| , p2[ for QCD 
SR. The latter measurement in [22] confirms the validity of the Gaussian ansatz for Ms(z 2 ) (at a 
small distance) as well as the value of the parameter A 2 . 

2 "Duality" transformation 

To obtain sum rule, we start with a 2-point correlator Uwxp (q) of tensor currents J^(x) = 
u(x)a^ v (zV) N d(x) {z is a light-like vector, z 2 = 0), 



n 



$?{Q) = *J d*x e^(0\T [J^(x)J^(0)\ |0) (4) 

whose properties were partially analyzed in ||, M. It is well known that the correlator at N = 
can be decomposed in invariant form factors n±, || ||| 

U^(q) = Il_(q 2 )P^ + U+(q 2 )P^ (5) 

where the projectors -Pi, 2, obeying the projector- type relations 

(Pi ■ P 3 Y^ = prvrprrv* = SijPf 1 ^ (no sum over i), P^ v = 3, (6) 

are presented in Appendix A. For the general case N ^ 0, a similar decomposition involves 4 new 
independent tensors Qi\ they appear due to a new vector z a introduced into the composite tensor 
current operator, 

U^ N f(q) = n-(q 2 ,qz)Pr aP +n+(q 2 ,qz)Pr aP +K 1 (q 2 ,qz)Qr af} 

+ K z {q 2 ,qz)Q^ + K z {q 2 ,qz)Qr aP + K q {q 2 iq z)Q^ . (7) 

Contributions of DAs, defined in Eqs.( [A.9"| )-( [A.10[ ), to different tensor structures in decomposition 
(0) are mixed, see Eqs.( A.ll )-( A. 12 ). The most effective way to disentangle them in practical OPE 
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calculations is to use explicit properties of different OPE terms under the duality transformation D 
(introduced in our previous work [16|) mapping any rank-4 tensor T^ u ' a/3 to another rank-4 tensor 

Dtf&p = X^'W* and & = L ( § ) 

Our projectors p^ u,a ^ ^ pV v -> a P Q^ u ' a/3 ^ Q% v , Q^ u ' al3 , and Q^ u '' a @ transform into each other under 
the action of D: 

(D Pl y^ = pr a " ; (dq^Y^ = [Pi + p 2 - Q 3 r a " oo 



DP, 



pr aP ; (D Q ^ aP = [p 1+ p 2 - Ql r^ ; do) 



(bQ z y u '^ = -Qr^; (bQ q y u '^ = [Q q -q z + q 1 + q 3 - Pi - p 2 \^ . (n) 



We have shown in |lq| that all terms in OPE could be divided into two classes, self-dual (DXgrj = 
Xgpj) and anti-self-dual (DX^sd = — ^ASd)- F° r example, the perturbative term is of ASD type, 
whereas the 4-quark scalar condensate contribution to OPE is of SD type. 

Below we introduce the shorthand notation for contributions of DAs to decomposition (|7j): vq, 
V\, and v 2 stand for 1 (p±, p' ± ); and uq, u\, and u 2 , for 1 H (bi), see Appendix A for details. For 
SD parts of OPE itj = — V{, and the system of equations simplifies to: 

n T (g 2 ,gz) K lt3 (q 2 ,qz) K q (q 2 ,qz)(=-2K z (q 2 ,qz)) 
2{qz)Nq 2 = Tvo -v t -v 2 ; 2{qz)Nq2 = T vi + v 2 ■ I( - jWq - 2 = v 2 , (12) 

whereas for ASD parts Ui = Vi, and we have: 

n T (g 2 ,^) K lt3 (q 2 ,qz) K z (q 2 ,qz) 2 

2{qz)Nq 2 = ™ + v 1+ v 2 ; 2{qz)Nq2 = - Vl - v 2 ; -^-^ = +v 2 ; K q (q ,qz)=0. (13) 

By these formulas, it is possible to determine p- and 6i -meson DA contributions of leading and higher 
twists. 



3 The "mixed parity" sum rule 

The usual way ||, |j to extract the moments of the function (f T (x) appeals to a correlator J(n,o)(q 2 ) 
of currents J^(0)z a and Jj^(x)z 13 defined as 

- M» (zq) N+ * J {Nfi) (q 2 ) = n{^(g) (V^) = "'^ ~ V ) 2 , (M) 

the latter equality in ( |14| ) follows from (|7|) and Eqs.( |A.7[ ) in Appendix A. This correlator contains the 
contributions from states with different parity, H-(q 2 ) and 11+ (g 2 ) (see the analysis in jl), therefore, 
the contamination from &i-meson in the phenomenological part of the corresponding 

SR is mandatory. The contamination makes it difficult to reliably extract the meson characteristics 
from this "mixed" SR. 

The main feature of the theoretical part of J(n,o)(q 2 ) is t ne cancellation of the self-dual part, 
represented by the four-quark condensate, in the anti-self-dual expression (|l~4|). The remaining "con- 
densate" parts of Eq. (O) contain, after the Borel transformation, the same 5 universal elements 
A<S> r (x;M 2 ) as for the p 7 S 7r-cases and, besides, an additional gluon contribution A<J>' G (x; M 2 ) (see 
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Appendix B). This term affects the values of moments rather strong, as was shown in ||. The con- 
tributions from the different kinds of NLC, A&y{x; M 2 ), are symbolically noted in the r.h.s. of SR 
(|i~5|). So, here we get rid of the four-quark condensate that is not known very well due to a possible 
vacuum dominance violation. But, the price we pay for it is a high sensitivity to an ill-known gluon 
contribution A& G (x; M 2 ). 

The method of calculation of the NLC contributions A<I>r {x\ M 2 ) to the theoretical part of SR 
is described in ^, [|]. The corrected final results of the calculation are presented in Appendix B 
that contains all the needed explicit expressions of A$r(^; M 2 ) for the simplest physically motivated 
Gaussian ansatz. The final SR including DAs of p-meson and next resonances p' and b\ into the 
phenomenological (left) part is as follows: 

2 IM 2 I b mivoJ / T 1 T\ _o/J1f2 



(/J) V(^" m?/M2 + (P "> P) + {fl)< i-)e- m ^' M = / " P f xed (», *; s T p , si) e^ M2 d s 

J 

+ A$ G (x; M 2 ) + A& G (x; M 2 ) + A$ v {x; M 2 ) + A<Z> T (x; M 2 ) , (15) 

where and are the effective continuum thresholds in p- and &i-channels. Recall again that 
the variation of the ill-known part of gluon contribution A& G (x; M 2 ) can reduce the second mo- 
ment significantly ||. In that paper, we suggest the following naive model: instead of the constant 
contribution Aip' G (x;M 2 ) = (a s GG) /(6nM 2 ) (as in the standard approach), we put 



A<f>' G (x; M 2 ) = Acp' G (x; M' 



0(A< x)6(x< 1- A) 



1 - 2A 

This simulation eliminates end-point (x = 0,1) effects due to the influence of the vacuum gluon 



nonlocality inspired by the analysis in [17] and our experience in the nonlocal quark case. The 
corresponding SR leads to estimate (£ 2 )p = 0.329(11) (see Fig. 2(a) ). However, this value drastically 

changes, (£ 2 ) T p -> 0.231(8), if we take the local expression Aip' G (x,M 2 ) unchanged. Therefore, the 

estimate (£ 2 )J = 0.329 contains a significant model uncertainty, and the real value seems to be 
smaller. 

Which prediction for this quantity can be obtained within the standard QCD SR approach? As 
one can see from Fig. 2(b) (long-dashed line), the value of (£ 2 )p cannot be estimated with a reasonable 
accuracy, because the standard SR does not have real stability. Nevertheless, the authors of Q 
bravely deduce an estimate (£ 2 )Jjb&B] = 0-27(4). We discuss this attempt in comparison with 
processing other SRs in greater detail in section 5. 



4 The "pure parity" sum rules 

Using the approach of Section 2, we calculate OPE terms for II^, K\^, and K z>q correlators and 
extract the contributions to DAs of the p- and &i-mesons. This allows us to write down the SRs for 
DAs of the p- and 6i-mesons separately: 

(m P fj)\ T p (x)e-^ M2 + (m p/ fJ,)\ T p ,(x)e- m l> /M2 = \ ff p^\x;s) s e'^ds 

+ A$ G (x; M 2 ) + A§s{x; M 2 ) + A<ly(x; M 2 ) + A<l T (x; M 2 ) ; (16) 

(rn b jl)\l{x)e-< IM2 = \ p^(x;s) s e'^ds 
+ A$ G (x; M 2 ) - A$ s (x; M 2 ) + A$ v (x; M 2 ) + A$ T (x; M 2 ) . (17) 
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where s^. b are the effective continuum thresholds in the p- and the 61-meson cases, respectively. The 

perturbative spectral density p^ rt (x;s) is presented in an order of 0(a s ) in |4], |8|] (Appendix B). 
Here we also define "tilded" functions 

A4> r (x; M 2 ) = ^M 4 <9 A/2 A$ r (*; M 2 ) , (18) 

and the whole tensor NLC contribution 

A<l r (x;M 2 ) = A$ Tl (x; M 2 ) + A$t 2 {x; M 2 ) - A$t 3 (x; M 2 ). (19) 

The later noticeably differs from the case of longitudinally polarized p-meson due to the opposite 
sign of T3-term, cf. ||. The theoretical "condensate" part in (|i~6|)-(p^|) contains 5 elements obtained 
from ( |i~8| ) with the same A$y{x;M 2 ) as for the p^-meson case, whereas the self-dual four-quark 
contribution A^s( x i M 2 ) is a new element of the analysis. Note, just this self-dual part A<l>s(x; M 2 ), 
entering in the SRs (|16|) and (0) with different sign, provides the different properties of the p- and 



&i-mesons fig] . 



For better understanding of the SR properties it is instructive to reduce them to standard version 
for (^ Ar )-moments. To this end, let us take the limits A 2 — > 0, A<3?r(x,M 2 ) — * Aipr(x, M 2 ) in 
eqs . (P^6j) - (P^?|) and integrate in x with weights (1 — 2x) N to obtain the local limit version of moment 
SR: 

(™,/j)V)Je-^ 2 = (20) 
1 [*J P e r t, . > (a s GG) f N - 1\ 16vr (^lg(O)g(O)) 2 



2 



i Orp (x:s) s e ' as — —77 

HT y ' 24vr \N + lJ 81 M 2 



S / M ds - - ) - w^si^m^j/ ^ N _ 



This SR demonstrates a considerably lower sensitivity to the gluon condensate contribution: the 
gluon part does not depend on the Borel parameter M 2 at all, and its relative value is 6 times as 
low as that in the "mixed" SR. The r.h.s. of Eq. ( |2(i| ) is reduced at N = to the known expression, 
see Q, that is not sensitive to the p' contribution, while its nonlocal version analyzed in |16| makes 
it possible to analyse the p' meson. For N > 0, the SR is unstable due to the effect of radiative 
corrections, and to obtain the moment estimates, we should return to the nonlocal version, Eq.([T6|). 

But the price one pays for this is high, the fidelity windows of the SRs are significantly reduced. 
For the /3-meson case, fidelity windows of the Borel parameters M 2 shrink to M 2 = 0.7 — 1.15 GeV 2 
(to be compared with M 2 = 0.75 — 2.25 GeV 2 in "mixed" SR) and demand one to take into account 
the //-meson explicitly. Here we cannot obtain the p'-meson mass from SR ( |l6|) because of the 
enhanced perturbative spectral density (~ s; this means that the differentiated SR has a spectral 
density ~ s 2 and presumably, is not stable at all); instead, we use the p'-meson mass extracted in 
our previous paper on the longitudinally polarized p- meson DA [P], m p i = 1496 ± 37 MeV, rather 



close to the Particle Data Group value m p i = 1465 ± 22 MeV []23 |. 

In the case of 61-meson, one can analyze only the SR for the zeroth moment (decay constant fjf ) 
of the DA (see Fig. 3), the SRs for higher moments appearing to be invalid. 



5 Processing different SRs and comparison of the results 

We start with considering the results of processing both the types of SRs for fj. Its dependence on 
the Borel parameter M 2 obtained from the "mixed parity" NLC SR, Eq. (|i~5|), with so = 2.9 GeV 2 
is shown in Fig. |l|(a). Figure ||(b) shows /J as a function of the Borel parameter M 2 obtained from 
the "pure parity" NLC SR, Eq. (H), with s = 2.8 GeV 2 . Both kinds of SRs are rather sensitive to 
the //-meson contribution and, for this reason, they were processed with taking it into account (see 
numerical results in Table 1). Solid lines correspond to the optimal thresholds sq] the dashed lines - 
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Figure 1: /J as a function of the Borel parameter M 2 obtained from: (a) the "mixed parity" NLC SR, Eq. with 
s = 2.9 GeV 2 ; (b) the "pure parity" NLC SR, Eq. with s = 2.8 GeV 2 . The fidelity windows for both figures 

coincide with the whole depicted range of M 2 . 



to the curves with the 10-fold variation of xLin (this corresponds approximately to the 5%-variation of 
so; definition of x 2 , see in Appendix C, Eq. QC.lp ). So, one can conclude that both types of NLC SRs 
agree rather well about the value of fT . Note that the presented fT is rather close to the standard 



estimation /J = 0.160(10) GeV |§ and to the lattice one /^ att (4GeV 2 ) = 0.165(11) GeV |2§], and 



differs significantly from the result /r = 0.140 GeV in p 







Table 1: The moments (£ N }m(p 2 ) at fi 2 ~ 1 GeV 2 
(errors are depicted in brackets in a standard manner) 


Type of SR 




f M (lGeV 2 ) 


N = 2 


N = 4 


N = 6 


N = 8 


Asympt. WF 


1 


0.2 


0.086 


0.047 


0.030 


NLC SR Eq.© : 


T 

P 


0.157(5) 


0.296(20) 


0.196(6) 


0.132(5) 


0.089(4) 


NLC SR Eq.© : 


T 

P 


0.162(5) 


0.329(11) 








B&B SR : 


P T 


0.160(10) 


0.304(40)[] 


does not work 


NLC SR Eq.© : 


P' T 


0.140(10) 


0.086(6) 


0.010(1) 


0.013(1) 


0.022(2) 


NLC SR Eq.© : 


hi 


0.184(5) 


does not work 


NLC SR Eq.© : 


h\ 


0.181(5) 


0.144(15) 








B&B SR : 


hi 


0.175(5) 


does not work 



Now we consider the results of processing SRs for the second moment (£ 2 )p . First, we demonstrate 
the results of the "standard" approach: (£ 2 )J from Eq.(3.21) in [Q] as a function of M 2 is shown in 
Fig. 2(b) by a long-dashed line. This curve is not stable in M 2 at all, therefore the SR can provide 

4 The estimate presented in this cell has been obtained by processing the "mixed parity" SR established in Q, 
whereas in the original paper fl] this value amounts to 0.27(4). 
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merely a range of admissible values, 0.27 < (£ 2 ) T p < 0.4. As it is evident from Fig.^, this wide window 
agrees reasonably with both the estimates from the "mixed" (a) and "pure" (b) NLC SRs. 



0.5 




0.75 



Y 

- \ 


(e) T P 


(b) 


\ \ 






\ 








— : 


0.6 


— 

0.8 1 1.2 1.4 1 



[GeV 2 ] 



Figure 2: (£ 2 ) T p as a function of the Borel parameter M 2 obtained from: 



M 2 [GeV 2 ] 

(a) the "mixed parity" NLC SR, Eq. (Q, 



with s = 2.9 GeV 2 ; (b) the "pure parity" NLC SR, Eq. @, with s = 2.8 GeV 2 . Both kinds of SRs were processed 
with taking the p'-meson into account, the arrows show the fidelity window (for figure (a) the window coincides with 
the whole depicted range of M 2 ). Solid lines correspond to the optimal thresholds so, the short-dashed lines on both 
the figures correspond to the curves with the 10%- variation of so (a) or of Xmin 3 )- The long-dashed line in figure (b) 
represents the SR of Ball-Braun M . 



Note, the authors of Q dealt with the quantity <22, the Gegenbauer coefficient in the expansion 
of DA. The second moment of DA is trivially connected with this coefficient, (£ 2 ) = 0.2 + (12/35)a2- 
Using the SR of Q for a<2, we obtain the corresponding window, 0.2 < 02 < 0.4, that leads to the 
mean value (£ 2 ) J pg tanc n = 0.30 being surprisingly close to our estimate from NLC SRs, (see Table 1). 



However Ball and Braun have obtained the erroneous estimate 02 = 0.2 ±0.1 producing, instead, the 

mean value (£ 2 >p [ B &B] = °' 27 - 

The curves for the next higher moments, whose estimates are presented in Table 1, have the 
fidelity windows and the stability behavior similar to (£ 2 )J (M 2 ) in Fig. 2(b). Finally, in Fig. 3, we 
demonstrate a very good correspondence between the values of obtained in different NLC SRs. 
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Figure 3: The curves in M 2 obtained from: (a) the "mixed parity" NLC SR (with taking the p'-meson into 
account with f p i defined from "pure parity" SR ©); (b) the ' 'pure parity" NLC SR (tljj). The arrows show the fidelity 
window (for the right figure, the window coincides with the whole depicted range of AP). Solid lines correspond to the 
optimal thresholds; the short-dashed lines on both the figures, - to the curves with the 10-fold variation of Xmini the 
long-dashed line on the right figure corresponds to the real B&B curve. 
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6 DA models and their check 



Possible models of DAs corresponding to the moments in Table 1 are of the form 



T.mod 

<Pp' 



(x,p 2 ) 



1.382 [p as {x)] z (l + 0.927Cf /2 (£) + 0.729Ct /z {0 



-,3/2, 



<p as (x) (1 + 0.29C 2 3/2 (e) + 0.41C7 4 o/z (0 



,3/2, 



0.32C 3/2 (£) 



T,mod/ 2 



(x) (l - 0.339C 2 3/2 (O + 0.003C 3/2 (£) + 0.192Cg /2 (£) 



<"W 2 ) = <p"{x)(l- (0.175 ±0.05)C 3/2 (O 



(21) 
(22) 
(23) 

where £ = 1 — 2x, C^(£) are the Gegenbauer polynomials (GP), and the norm p 2 ~ 1 GeV 2 corre- 
sponds to a mean value of M 2 . Recall again that the value of the important coefficient a2 = 0.29 in 
(21) is confirmed by 3 sources: "pure" NLC SR ([U|), "mixed" NLC SR flUD, and a mean value from 
the "mixed" standard SR. Figures ||(a) contain curves of DA corresponding to p±, eqs. (|2l|), and 
p' ± (22). The arising 3- hump shape of DA for p± drastically differs from that obtained in Q and 
from the one obtained in chiral effective theory . 




Figure 4: (a) The curves of ipJ' mo (a;, 1 GeV 2 ): Solid lines correspond to the best fits for determined moments (see 



Table Vh the dashed line on the left figure corresponds to the B&B curve (which fits only {£ 



in x. Different lines here correspond to different values of Borel M 2 — 0.7 — 0.9 GeV 



0.27). (b) The rhs 





X 



X 



Figure 5: (a) The curve of ip^; mod (x, 1 GeV 2 ) in x. (b) The rhs of Eq.@ SRj(a;,M 2 ) in x. Solid and dashed lines 
here correspond to different values of nonlocality parameter A 2 = 0.4 — 0.5 GeV 2 with fixed value of Borel parameter 



0.8 GeV 



This difference mainly appears due to the higher moments, A = 4, 6, 8, involved into considera- 
tion. Nevertheless, the hump shape is not an artifact of the GP expansion series truncation. These 
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models really contain only 3 first GPs, meanwhile, it is enough to reproduce all 4 moments up to 
N = 8. Moreover, an additionally smoothed^ rhs of the NLC SR (|l6|) demonstrates qualitatively 
the same behaviour in x (at admissible M 2 ) as the model DA, compare Figs. |](a) and (b). The 
stability of the DA shape with respect to the variation of ansatz is also checked. To this end, we 
show in Fig. ||(b) the same r.h.s. of ( |i~6|) as in Fig. |](b), but with different values of the single ansatz 
parameter A 2 = 0.4 - 0.5 GeV 2 at fixed value M 2 = 0.8 GeV 2 . 

Inverse moments of DAs often appear in perturbative QCD predictions for exclusive reactions. 
The estimates for important (x~ 1 )m moments obtained from the model DAs are presented heref] 





[x, 1 GeV 2 ) 




x 




[x, 1 GeV 2 ) 




X 




[x, 1 GeV 2 ) 


X 



dx 

dx 
dx 



4.15 



-0.4 
-0.1 



[ 3.6 
2.57 ± 0.20 (here 

2.48 ± 0.20 (here 



(here) 

(B&B model) 



(24) 

(25) 
(26) 



It is useful to construct an independent SR for these inverse moments to verify the DA models 
( |2~T| . [22] , |23| ). Namely, the weighted sum C(M 2 ) of these moments 




C(M 2 ) = {x^) p + {x^) p , ( J JL ) e -K'-^)^ 2 + <,.-),„ 
can be obtained by integrating the rhs of the "mixed" NLC SR ( |15[) with the weight 1/x 




(ml -m 2 p )/M 2 



(27) 




. 8 



1.2 1.4 1 . 



1 . 1 



Figure 6: C(M 2 ) as a function of M 2 (solid line) determined by Eq.(|27| and integrating in x Eq.(^H|) in comparison 
with the lhs of Eq.(pi|) (long-dashed line). The dotted line corresponds to (x~ 1 ) p — 4.15; whereas the dashed line - to 
the lhs of Eq.(^) with upper values of corresponding moments. 

A comparison of the function C(M 2 ) with the corresponding combination of model estimates 
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) obtained in different kinds of NLC SRs (mainly from the "pure" ones) leads to an approximate 
equation 



4.15 + 2.57 




e - {m l'~ m2 " )/M2 +2A8 




"K, -m 2 p )/M 2 



C(M 2 



(28) 



illustrated in Fig. ^. 

As a result, one can conclude: 



5 A certain smoothing of some (^-functions in the r.h.s. of the SR (see Appendix B) is not important. 
6 The upper error +0.4 in (E4J) corresponds to an overestimate (£ 2 ) = 0.329 from the "mixed" SR 
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1. The "mixed" NLC SR is highly sensitive to b\- and p'-meson contributions, the difference in 
the behavior of C(M 2 ) (solid line) and in the /^-contribution alone (dotted line) illustrates this 
point. 

2. The curve C(M 2 ) lies between mean and upper estimates for the lhs of d28|), so it is in reasonable 
agreement with the estimates (0, |2 

moments in the "mixed" SR as compared to that obtained from the "pure" one. 



26). It also demonstrates an overestimation of DA 



7 DA models and the B —> peu decay form factors 

The new DA shapes result in different pQCD predictions for exclusive reactions with the p-meson. As 
an example, we re-estimate form factors V(t), A\^{t) corresponding to the transition matrix element 



(p,X\(V — A)a\B) of the process B — > p eu, in the framework of the light-cone SR approach, [26|. 
That was done earlier by Ball and Braun in p7|| , [ p£| ] on the base of DAs from Thus, to estimate 
the influence of the new nonperturbative input presented in the previous sections, we have used the 
LC SR in the leading twist approximation (cf. |27]]). Just as in the case of the LC expansion for the 
transition amplitude 7*7 — > 7r°, one might expect high sensitivity to the end-point behavior of the 



DAs, as they enter into convolution integrals like {x )m estimated in (24). 

However, there are some essential differences which effectively soften our expectations. First, the 
DAs also enter into the "phenomenological" side of the SR in the "continuum" contribution of higher 
excited states in the channel with l?-meson quantum numbers. This, actually, is a specific feature 
of any LC SR. By subtracting the "continuum" , one actually obtains "infrared safe quantities" like 
jl dxip(x)/x where e ~ (to 2 — t)/(s B — t), nib — 4.8 GeV, and s B ~ 34 GeV 2 is the continuum 
threshold in the .B-channef] as defined from the 2-point QCD SRs for the £>-meson decay constant 
/b (see |2^|). For t ~ 0, e ~ 0.5 — 0.6 and the LC SR should not be so sensitive to the end-point region 
x ~ 0. Obviously, the end-point region becomes to be important for higher momentum transfers t. 
However, for t > 20 GeV 2 the LC expansion would hardly make sense. The second factor which 
eventually decreases the importance of the end-point region is connected with the standard Borel 
transformation of the SR with respect to the virtuality of the 5-meson current: -p% -> M|. The 
corresponding contribution from the coefficient function produces a standard suppression exponent: 
exp(x(t — mfy/xMg). Numerically, it occurred to be less important. 

We have treated the LC SRs using the same input parameters and the same procedure of extract- 
ing the physical form factors as in Ref.[^]. However, if one tries to fix the onset of the "continuum" 
by hand to the value s B ~ 34 GeV 2 dictated by the 2-point SRs for one encounters inadmissi- 
ble uncertainties in the determination of the form factors when using our new nonperturbative input 
DAs. To get a stable SR, one is forced to allow a higher value for the s B parameter. 

Below, the form factor values are written at a zero momentum transfer (t = 0) as compared with 
B&B results: 

V(0) = 

^l(0) = \ n 07mV VlFVTtl rL U - Q/l n^2l J ' i ^ 1 1\ ( 29 ) 



0.37(1) 


(here 


[so = 


50 GeV 2 ], 


x 2 - 


i 0.4) 


0.35(2) 


([27] 


[so = 


34 GeV 2 ], 




3.4) 


0.283(4) 


(here 


[so -- 


= 45 GeV 2 


, x 2 


« 0.1 


0.27(1) 


(H 


[so = 


= 34 GeV 2 ] 


, x 2 


* Lb 


0.30(1) 


(here 


[so = 


50 GeV 2 ], 


x 2 * 


i 0.2) 


0.28(1) 


(H 


[so = 


34 GeV 2 ], 


x 2 * 


1.1) 



^4 2 (0) = 

Our form factors are slightly higher than those in |^7j and possess a better accuracy (compare \ 2 
in (|30|)). The difference becomes more pronounced for a large value of the momentum transfer t, 

7 As we shall see below, the LC SRs "prefer" a higher value. 
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(m? — t ~ 0(mb)). The last is not surprising due to higher sensitivity to the end-point behavior of the 
input DA in this region. The form factors presented are determined with new "optimal" thresholds 
Sq providing few times better processing accuracy. Note that the parameters of the usual "pole" 
parameterization of the form factors change significantly as compared to that in [27], e.g., 



. , , 0.283 



1 - 0.157 (t/m%) - 0.837(t/m|) 2 

The important form factor A\(t) (solid line) increases about 5 — 10% in comparison with the B&B 
result (the bars in the figure show the errors of the B&B calculations), with an optimal threshold 
s£ ~ 45 GeV 2 . 




Figure 7: Form factor Ai(t); solid line corresponds to our processing of LC QCD SR, dashed line - to processing 
following the B&B formulas (ETA) (the bars in the figure show the errors of the B&B calculations). 



8 Conclusion 

Let us summarize the main results of this paper: 

1. We construct NLC SRs for DA for each P-parity channels, based on the properties of the duality 
transformation. The negative parity NLC SR for transversely polarized p-, p'-mesons works 
rather well and allows us to estimate the 2-nd, 4-th, 6-th, and 8-th moments of the leading 
twist DAs. The positive parity SR for the transversely polarized 61-meson can provide only the 
value of the 61-meson lepton decay constant, . It should be emphasized that an analogous 
evaluation of the moments within the standard QCD SR approach is impossible. 

2. Results of processing different NLC SRs of the "pure" (see Figs, lb, 2b, 3b) and "mixed" (see 
Figs, la, 2a, 3a) parity are compared, and a reasonable agreement between them is found. The 
"mixed" SR in the standard version admits merely a window of possible values of the second 
moment (£ 2 ) (see, e.g. , [Q]); the position of the window is corrected here and, as a result, agrees 
with the NLC SR results presented in Table 1. 

3. The models for the leading twist DAs of the p±- and //j_- m esons, ( p|j22] ), and of the fr^-meson, 
(HH), are suggested. The shape of a new /?j_-meson distribution (see Fig. 4a) drastically differs 
from that obtained by Ball and Braun Q only on the basis of the value 02 = 0.2. The latter 
estimate is discussed in sect. 5. 
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4. We estimate important integrals appearing in perturbative QCD predictions for different ex- 

elusive reactions, {x~ 1 )m = / dx in (24)-(p3), based on our results for the DA shapes. 

Jo x 

We check the self-consistency of these results by comparing them with those obtained from an 
independent "mixed" QCD SR for the inverse moment {x~ 1 )m and find an agreement. 

5. Form factors of the process B — > p ev, V(t), Ai^(t) where t is momentum transfer are also 



re-estimated in the framework of the light-cone SR approach [27| on the basis of the new model 
for the p-meson DAs; the results are slightly higher and have uncertainties a few times as small 
as those obtained by Ball and Braun. 

Finally, we can conclude that the nonlocal condensate QCD SR approach to distribution ampli- 
tudes is self-consistent and gives reliable results. An open problem of this approach is to determine 
well-established models of distribution functions /r(^) from the theory of nonperturbative QCD vac- 
uum. First direct attempts to calculate quark NLC have been done in lattice simulations in pT| . 
The "short distance" correlation length of NLC has also been extracted later in [22|; it turns out to 



be reasonably close to the value of 1/X q and confirms the validity of our Gaussian NLC model. 
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Appendix 

A Decomposition of rank-4 tensor II^ a/3 



Q/xu;a/3 



1 

v 



gi ia q u ql 3 — g ua q^q( } — g^q u q a + g v ^q^q a 



fiv;af3 



1 



2(qz) 
1 

„2 



gTtfzP + g vfi q^z a - g^q v z a - g^q^z? 



2{qz) 2 



Z$ - qPz a 



(A.l) 
(A.2) 
(A.3) 

(A.4) 
(A.5) 
(A.6) 



g^ a 'z l ' z^ 'Pf y,a ^ = p!f z '^ z = —p^ z '^ z 



S 2 



; QT' 



qfq a z u z (3 Pr ,aP = PT' qZ = QT qZ = QT qZ = -Qf' qz 



3 Q 
(qzf 



fiz;/iz 



. r>qz;qz 
i r 2 



= Q9m» = o. (A.8) 
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Let us write down the parameterization of matrix elements of a composite tensor current operator, 
see, e.g. , @: 



(0 | d(z)(v«(0) | p±(p,A)> 



z 2 =0 



(eJp, X)p u - e u (p, X) Pli ) / dx <p£{x) e ix( - zp) 



o 



+ (e„(p, A>„ - e„(p, A)^) p 2 / Vi(x) e*^ 



+ - p%) (e(p, X)z)p 2 / dx V 2 (x) e ix ^ 



(A.9) 



(0 | d{z)a^u(0) | 6i(p,A)) 



2 2 =0 



W £ >> A W dx <p h (x) e ix ^ 



ix(zp) 



+e^vap£ a (p,X)zPp 2 I dx Ui(x) e 



-e lluaf ip a zf 3 (e(p,X)z)p 2 / dx U 2 (x) e**™ 



(A.10) 



Here we decode our shorthand notation used in Section 2 

2 



u 



fp± ( xN )p± ; v ± 



fj ± 2 {-iNx N ~ l ) Vl ■ v 2 = fj ± 2 (-N(N - l)x N ~ 2 ) V2 ; 



fb ± 2 (x N )b x ; «i = £ 2 {-iNx N ~ 1 ) Ul ; u 2 = /£ 2 (—N(N - l)x N ~ 2 ) U2 



(with {f{x))u = Jq dx f(x)U(x)). In the general case, the whole system of equations for different 
twist DA contributions is of the following form 

n-(q 2 ,qz) , , ifi(g 2 ,gz) K z {q 2 ,qz) 

- -v + ui + u 2 ; 7^—^r5. = -vi - u 2 ; n/ _._ AN _. 9 = +u 2 ; (A.ll) 



2(qz) N q 2 
n + (g 2 ,gz) 



2(g*)V 



K 3 (q 2 ,qz) , , ,,,, 

+UO + U1+U2; „^... A ^.. 9 = -^1 - U 2 ; TTTTTwvTcT = u 2 ~ U 2 . (A.12) 



2(^)^2 
Kq(q 2 ,qz) 



2(qz) N q 2 - , ^ . ^ , 2(^) V " - - > 2(^)^2 

B Expressions for nonlocal contributions to SR 

To construct SR for distribution amplitudes, it is useful to parameterize NLC behaviors by the 
"distribution functions" ||, [14| a'la a-representation of propagators, e.g. , fs(a) for the scalar 
condensate Ms{z 2 ) | 

M S (z 2 ) = (g(O)g(O)) jf e^ 2 / 4 f s (a) da, where jf /<?(«) da = 1, J af s (a)da = -j-, (B.l) 

and for the vector condensate My(z 2 ), 

MMz) = (q(0)^q(z)) = -iz» / e az /4 f v (a) da, where / / v (a)da = 0. (B.2) 

4 Jo J 

Here and in the following we take quark and gluon fields in the fixed point gauge z^A^(z) = 
where the path-ordered exponential E(0, z) = 1. The appearing in the SR quark-gluon-quark NLC 

M^(M^), 



M^(y, z) = (q{<l) lu A^z)q{y)) = (y^z v - g^zy)) ■ M Tl + (z^ - g^z 2 ) ■ M T2 + ■ • <B.3) 
M^iy, z) = (q{0hv{7 5 )M z )q(y)) = ^upaZ p ya • M T3 + • • • , (B.4) 



In deriving these sum rules we can always make a Wick rotation, i.e., we assume that all coordinates are Euclidean, 



z z < 0. 
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can be decomposed in form factors Mti-T3, where the tensors in front of them satisfy the gauge 
condition z^M^ LV (M lJliU ) = ( since z^A^(z) = 0). The NLC Mti-tz can be parameterized by a 
triple integral representation 



Mn(z ,y , (z — y) ) = A 



; ( ttlz V4+<» 2 s J /4+«3N 8 /4)j.( Ql)a2i a3 ) da 1( la 2 da 3 , 



(B.5) 



8tt 



where At% = { — §As> \A$, § As}) and As = — (^/a^q(0)q(0)} 2 . The function /5(a) and other similar 
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functions /r(«) describe distributions of vacuum fields in virtuality a for every type (T) of NLC. 
The convolutions A&r(x, M 2 ) of the distribution functions fr and coefficient functions completely 
determine the r.h.s. of SR's, so A$r depends on the model of fr- For vacuum distribution functions 
fr(oi), we use the set of the simplest ansatzes 



fs(a) = s(a-X 2 q /2) ; f v (a) = 5'(a-\ 2 q /2) ■ 
/ti, 2 , 3 («i 5 «2,Q!3) = 5 (ai - Xg/2j 5 (a 2 - X 2 q /2j 5 (a 3 - \ 2 q /2 



(B.6) 
(B.7) 

Their meaning and relation to initial NLCs have been discussed in detail in |], f?J. The contributions 
A$t{x, M 2 ) to the r.h.s of SR, corresponding to these ansatzes, are shown below. The limit of these 
expressions to the standard (local) contributions fr(x, M 2 ) - X 2 — > 0, A&r(x, M 2 ) — > Apr(x, M 2 ) 
are also written for comparison. Hereafter A = X 2 /(2M 2 ), A = 1 — A: 

A s 18 



A$ 5 (x,M' 



A<ps 
A$ v 
Aip v 
A$ Tl 

Aip Tl 
A$t 2 

A(p T2 

A$t 3 

Aip T3 

A&g 



x,M 2 
x,M 2 
x,M 2 
x,M 2 

x,M 2 
x,M 2 

x,M 2 
2 



M 4 AA 2 



A 



{6 (x > A > x) x [x + (A — x) In (a;)] + (x — * x) + 
+9(1 > A)9 (A > x > A) [A + (A - 2xx) ln(A)] } , 



^9(<5(x) + (x^x)); 
^(x5'(x-A) + (x^x)), 
-j-fi (xS 1 (x) + (x -> x)) ; 

M 4 



;i > 2A) j[5(x - 2A) - 5(x - A)] ~ 2 ) + d ( 2A 



> xY 



x 



6(x>A)^ 



x-2A 



AA 



+ (x 



(B.8) 

(B.9) 
(B.10) 

(B.ll) 



3A S 

4A S 



(6' (x) + (x — > x)) ; 



X0( 



M 4 
1 + 2x - 4A 



1 > 2A){ ^ X A 2A) - #(2A > x)6(x > A)- 



+ (x 



A 



2A S 



M 4 



AA 2 

(x<5' (x) + (x — > x)) 



(B.12) 



= T^nr l»(2A > > A)0(1 > 2A) 



x,ikP 



x,ikP 



M 4 AA I 

+ (x — > x) 
3A S 



x A 

2 -A-A 



M 4 

(a s GG) 



2AnM 2 



(<5(x) + (x^x)); 

(S(x - A) + (x-> x)) 



(B.13) 



(B.14) 
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Aip G (x, M 2 
A& G (x,M 2 
Aip' G (x, M 2 



(a s GG) 



(5 (x) + (x — > x)) ; 



24vrM 2 

(a s GG) d(A < x)6(x < 1 - A) 



6vrM 2 
(a s GG) 



1 - 2A 



(B.15) 



6vrM 2 " 

For quark and gluon condensates, we use the standard estimates {y/a^q(0)q(0)) ~ (—0.238 GeV) 3 , 

= 0.4 ±0.1 GeV 2 normalized at \i 2 ps 1 GeV 2 . 



(asGG) - 0.001 GeV* @ and A 2 - G^ ** 



12tt ^ q 2( 

Expressions for perturbative spectral density: Radiative corrections reach 10 % of the 
Born result at s ~ 1 GeV 2 . 



4tt 



2 In 



7T 



+ 6 — + ln 2 (x/x) + ln(xx) 



(B.16) 



Here y? ~ 1 GeV 2 corresponds to the average value of the Borel parameter M 2 in the stability 
window; a, (lGeV 2 ) « 0.52. We also use the 'mixed' perturbative spectral density suggested in |3l| 
in the "mixed" SR: 



„mixed 

Pt 



(x, s; sj, si) = p^ ert (x; a)- [o (s T p - s) + 6 (sj - s 



(B.17) 



C About x 2 -definition in Sum Rules 

Let us discuss the definition of x 2 f° r the SR case. We have here the function F(M 2 ,s), and the 
problem is to find the best value sq, such that F(M 2 ,sq) is the most close to a constant value for 
M 2 < M 2 < (values of are known and fixed from standard constraints of QCD SR, see 



). We define the function x 2 ( s ) f° r tne curve F(M 2 , s) with M 2 E [M 2 ,M 2 ] in the following 
manner: 



[30 



X 2 (s) 



1 



(N - l)e 



N 

E 

k=0 



1 N 

F (Ml + kS,s) - —— J2f(mI + k5,s 

+ fc=0 



(C.l) 



where 5 = (Ml — Ml)/N, N ~ 10, and e is of an order of the last decimal digit in F(M 2 , s) we are 
interested in (in the case of decay constant f p ?s 200 MeV, e ?s 1 MeV; in the case of the second 
moment (£ 2 ) p w 0.25, e ?s 0.01). Then, if we obtain x 2 ( s o) ~ 1? this tells us that the mean deviation 
of F(M 2 ,sq) from a constant value in the region [Ml, Ml] is about e. To find the minimum value 
of x 2 and the corresponding sq, we used the code Mathematica. 



References 

[1] V. L. Chernyak and A. R. Zhitnitsky, JETP Lett. 25 (1977) 510; 

V. L. Chernyak, A. R. Zhitnitsky, and V. G. Serbo, JETP Lett. 26 (1977) 594; 
Stanley J. Brodsky and G. Peter Lepage, Phys. Lett. B87 (1979) 359; 
A. V. Efremov and A. V. Radyushkin, Phys. Lett. B94 (1980) 245. 

[2] V. L. Chernyak and A. R. Zhitnitsky, Nucl. Phys. B201 (1982) 492; B214 (1983) 547(E); 
Phys. Rept. 112 (1984) 173. 

[3] J. Govaerts, L. J. Reinders, F. de Viron, and J. Weyers, Nucl. Phys. B283 (1987) 706. 

[4] Patricia Ball and V. M. Braun, Phys. Rev. D54 (1996) 2182. 



16 



[5] Patricia Ball, V. M. Braun, Y. Koike, and K. Tanaka, Nucl. Phys. B529 (1998) 323. 

[6] S. V. Mikhailov and A. V. Radyushkin, JETP Lett. 43 (1986) 712; Sov. J. Nucl. Phys. 49 (1989) 
494. 

S. V. Mikhailov and A. V. Radyushkin, Phys. Rev. D45 (1992) 1754; 
A. P. Bakulev, A. V. Radyushkin, Phys. Lett. B271 (1991) 223. 

[8] A. P. Bakulev and S. V. Mikhailov, Phys. Lett. B436 (1998) 351. 

[9] A. V. Radyushkin, Talk given at "ICTP Conference on Perspectives in Hadronic Physics, Tri- 
este, Italy, 12-16 May, 1997" [|hep-ph/9707335| , pp. 1-10. 



[10] V. M. Braun and I. E. Filyanov, Z. Phys. C44 (1989) 157. 

[11] N. G. Stefanis, W. Schroers, and H. C. Kim, Phys. Lett. B449 (1999) 299. 

[12] A. Schmedding and O. Yakovlev, Phys. Rev. D62 (2000) 116002. 

[13] A. V. Radyushkin, in Workshop on Continuous Advances in QCD, Proc. Workshop, Minneapolis, 
Feb 18-20, 1994 (University of Minnesota, Minneapolis, 1994), pp. 238-248. 

[14] A. P. Bakulev, S. V. Mikhailov, Z. Phys. C68 (1995) 451; Mod. Phys. Lett. All (1996) 1611. 

[15] A. V. Radyushkin and R. T. Ruskov, Nucl. Phys. B481 (1996) 625. 

[16] A. P. Bakulev and S. V. Mikhailov, Eur. Phys. J. C17 (2000) 129. 

[17] S. V. Mikhailov, Phys. Atom. Nucl. 56 (1993) 650. 

[18] A. G. Grozin, Int. J. Mod. Phys. A10 (1995) 3497. 

[19] V. M. Belyaev, B. L. Ioffe, Sov. Phys. JETP 57 (1983) 716; 

A. A. Ovchinnikov, A. A. Pivovarov, Sov. J. Nucl. Phys. 48 (1988) 721; 
A. A. Pivovarov, Bull. Lebedev Phys. Inst. 5 (1991) 1. 

[20] M. V. Polyakov, C. Weiss, Phys. Lett. B387 (1996) 841; 

A. E. Dorokhov, S. V. Esaibegian, and S. V. Mikhailov, Phys. Rev. D56 (1997) 4062. 

[21] M. D'Elia, A. Di Giacomo, and E. Meggiolaro, Phys. Rev. D59 (1999) 054503; 
E. Meggiolaro, Nucl. Phys. Proc. Suppl. 83 (2000) 512. 

[22] H. G. Dosch, M. Eidemiiller, M. Jamin, and E. Meggiolaro, JHEP 07 (2000) 023. 

[23] D. E. Groom et al, Eur. Phys. J. C15 (2000) 1. 



[24] D. Becirevic et al, Nucl. Phys. Proc. Suppl. 73 (1999) 222; |iep-lat/9809129H . 
[25] M. V. Polyakov, Nucl. Phys. B555 (1999) 231. 

[26] A. P. Bakulev, S. V. Mikhailov, and R. Ruskov, in Heavy Quark Physics, Proc. 5th Int. Work- 
shop, Dubna, Russia, 5-9 Apr 2000, edited by M. A. Ivanov, V. E. Lyubovitskij, and E. Lipartia 
(Joint Institute for Nuclear Research, Dubna, 2000), pp. 31-35. 

[27] Patricia Ball and V. M. Braun, Phys. Rev. D55 (1997) 5561. 

[28] Patricia Ball and V. M. Braun, Phys. Rev. D58 (1998) 094016. 



17 



[29] T. M. Aliev and V. L. Eletsky, Sov. J. Nucl. Phys. 38 (1983) 936; 

E. Bagan, Patricia Ball, V. M. Braun, and H. G. Dosch, Phys. Lett. B278 (1992) 457. 

[30] M. A. Shifman, A. I. Vainshtein, and V. I. Zakharov, Nucl. Phys. B147 (1979) 385, 448, 519. 

[31] V. M. Belyaev and A. Oganesian, Phys. Lett. B395 (1997) 307. 



18 



